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ABSTRACT 

Wc revisit the anisotropic universe model previously developed by Ackerman, Carroll and Wise 
( ACW) , and generalize both the theoretical and computational framework to include polarization and 
various forms of systematic effects. Wc apply our new tools to simulated WMAP data in order to 
understand the potential impact of asymmetric beams, noise mis-estimation and potential Zodiacal 
light emission. We find that neither has any significant impact on the results. We next show that the 
pr eviously reported ACW signal is also present in the 1-year WMAP temperature sky map presented 
bv lLiu fc Lil (|2009af ). where data cuts are more aggressive. Finally, we reanalyze the 5-year WMAP 

data taking into account a previously neglected -term in the signal covariance matrix. We 

still find a strong detection of a preferred direction in the temperature map. Including multipoles 
up to i — 400, the anisotropy amplitude for the W-band is found to be g = 0.29 ± 0.031, nonzero 
at 9cr. However, the corresponding preferred direction is also shifted very clos e to the ecliptic poles 
at (7,6) = (96,30), in agreement with the analysis of lHanson fe Lewis! (|2009f) . indicating that the 
signal is aligned along the plane of the solar system. This strongly suggests that the signal is not of 
cosmological origin, but most likely is a product of an unknown systematic effect. Determining the 
nature of the systematic effect is of vital importance, as it might affect other cosmological conclusions 
from the WMAP experiment. Finally, we provide a forecast for the Planck experiment including 
polarization. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

In recent years, the study of the cosmic microwave 
background (CMB) has proved to be the most fruitful 
addition to our understanding of the early universe. Ob- 
servations of the CMB anisotropies, like those obtained 
by the Wilki nson Microwave Anis o tropy Probe (WMAP ) 
experiment ((Bennett et all 120031 : iHinshaw et all I2007D . 
have provided us with incomparable insight on the com- 
position of structure in our universe. Combined with 
previous experimental knowledge and a sound theoret- 
ical framework, the concordance model of ACDM has 
been established. 

The ACDM model relies on the framework of in- 
flation. Inflation was initially proposed as a solu- 
tion to the horizon and flatness problem (jGuth et al I 
Il981f ). Additionally, it established a highly suc- 
cessful theory for the formation of primordial den- 
sity perturbations, providing the required seeds for 
the large-scale structures (LSS). Eventually, these 
later gave rise to the temperature anisotropies in 
the cosmic microwave background radiation that we 
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observe todav (IGuth et al 


I198U iLinde et all 
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11982: 
PJ82; 


Mukhanov & Chibisov 11981 


Linde et all 19831. 1994 iSmoot et al. 


19921: IRuhl et al 


20031 iRvnvan et al.ll2003t IScott et all 


20031). 



One of the predictions from inflation is that the ob- 
served universe should be nearly isotropic on large scales. 
Howev er, anomalies found in the CMB during the recent 
years tide Oliveira-Costa et al.l 120041 : iVielva et all 120041 : 
lEriksen et al.ll2004al) suggest that anisotropic inflation- 
ary models should be considered. A specific example 
is the generalized model presented by I Ackerman et al.l 
(|2007t ). which considers violation of rotational invariance 
in the early universe. A general framework for describing 
simila r models was presented bv lPuUen fc Kamionkowskil 

(120071). 

iHimmetoglu et alj (|2009al lrJ) showed that the 
anisotropic inflationary background of the ACW 
model characterized by a fixed-norm vector field ulti- 
mately is unstable. However, the paramctrization of 
the signal covariance matrix is independent of that 
and is very useful. It represents general 
by rotations in the CMB at a 
Several papers have recently 
of the ACW model with 
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signal 
unstable model, 
correlations induced 
phcnomcnological level, 
investigated the properties 



extensions dHou et alJ 20091: iKarciauskas et al.l 120091 : 



iDimopoulos et al.ll2009t ICarroll et all 120081 ). 

Work in this field suggests that the 5-year WMAP 
data contains a significant ACW anisotropic signal, 
corresponding to a 3.8a de tection in the W-band 
(iGroeneboom fc Eriksenl 120091). A more recent paper by 
lHanson fc Lewis! (|2009P ) points out that the direction is 
incorrect due to a n eglected factor o f (— corrected 
in a later version of I Ackerman et all (|2007T ). yielding an 
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ACW-signal in which the preferred direction is located 
very close to the ecliptic poles. 

In this paper, we re-analyze the 5-year WMAP data 
including the previously neglected -factor, and 

investigate whether traces of the ACW anisotropic con- 
tribution signal are still evident. The analysis will, as 
previously, be p erformed with the CMB Gibbs sam- 
pling framework ll Jewell et al.l l2q04l;lWandelt et al.ll2004l: 
lEriksen et al.ll2004bh . which bv TGroeneboom fe Eriksenl 
(|2009f ) was included to allow for non-diagonal, but sparse 
covariance matrices. This framework allows for exact 
Bayesian analysis of high-resolution CMB data with 
a non-diagonal CMB signal covariance matrix. The 
isotropic method has already been applied several times 
to the WMAP data (|Q'Dwver et al.ll2004l : lEriksen et all 
l2007alTbl . I2008bf ) , and has already been exten ded to take 
into account polarization llLarson et al.ll2007l) a nd inter- 
nal component separation (jEriksen et al .I l2008alh 

In our re-analysis of the 5-year WMAP temperature 
data, we confirm that the direction is shifted to the eclip- 
tic poles, at a greatly increased significance. As the north 
and south ecliptic poles are aligned with our solar system, 
the ACW signal in the WMAP data is therefore most 
likely a systematic effect and not of cosmological ori- 
;in. T his is in complete agreement with lHanson fc Lewis! 
2009h . 

It has not yet been possible to fully rule out whether 
any known systematic effect could have contributed to 
the signal. In theory, either asymmetric beams, mis- 
estimated noise or even the Zodiacal light could have 
affected the detection of the ACW-signal. In this paper, 
we consider these three effects and conclude that neither 
have any effect on the ACW-signal. 

Until now, the framework has only supported 
temperature-temperature correlations (TT). Here, we ex- 
tend the mechanics to include E-mode correlations (EE), 
including cross- mode correlations (TE) . We then provide 
a forecast for the upcoming Planck experiment, consid- 
ering simulated T+E maps. The Planck data will hope- 
fully be able to rule out all doubts about the origin of 
the ACW signal. 

2. THE ACW MODEL WITH POLARIZATION 

We are interested in the signatures that the ACW 
model would leave on the polarization of the CMB and 
focus our attention on the scala r perturbations. This cal- 
culati on was first performed bv lPullen fc Kamionkowskil 
(|2007| ). Observing the CMB sky in the direction e pro- 
vides information of the E-mode polarization constructed 
from the Stokes parameters Q(e) and U (e), as well as the 
temperature T(e). One can express the respective maps 
in terms of the spherical-harmonic coefficients a,E,i m and 
a>T,im which are given for each X = {E, T} by 



mode and satisfies 

(<5(k) ( 5*(q))=P'(k) ( 5 3 (k-q). 



(2) 



The ACW model proposes that if we drop the assump- 
tion of statistical isotropy by having a preferred direc- 
tion n during inflation, the primordial power spectrum 
at leading order has the form 



P'(k) =P(k) (l+5(fc)(k-n) : 



(•3) 

Here, g(k) is a general function of fc, which ACW argue 
is well approximated by a constant, g*. 

To study the statistics of the CMB produced by the 
scalar perturbations we need the power spectrum of the 
T, E modes and the cross-correlation between them. Us- 
ing the expressions JT]) and ([3]) we can write the various 
correlations for X = {E, T} as 

{ax,lm a *X',l'm') = Sll'S mm i C^ X +g*£lm;l<m< Cfc? , (4) 

where the C^, x are given by 

/>oo 

H) 1 -'' / dfcfc 2 p(fc)e£)(fc)e^,(fc). (5) 



c, 



XX' 



I.I' 



The coefficients £, 



lm;l'm' encode the departure from 
isotropy and connect / wi th I' = {1,1 ± 2| and ro with 
m' = {m, m ± 1. m ± 2} (jAckerman et al.l l2007h . Note 



that the factor of (—i 
of the paper. 



was missing in the first version 



3. THE POLARIZED ANISOTROPIC CMB GIBBS SAMPLER 

CMB data observations can be modeled as: 

d = As + n, (6) 

where d represents the observed data, A denotes 
convolution by an instrumental beam, s(#, <p) = 
Y^e m a imXimiQ i </>) is the CMB sky signal represented in 
cither harmonic or real space and n is instrumental noise. 
It is generally a good approximation to assume both the 
CMB and noise to be zero mean Gaussian distributed 
variates, with covariance matrices S and N, respectively. 
In harmonic space, the signal covariance matrix is de- 
fined by S(.m,t'm' = {aimO*ei m ') ■ in the isotropic case, 
this matrix is diagonal. The connection to cosmological 
parameters u> is made through this covariance matrix. 
Finally, for experiments such as WMAP, the noise is of- 
ten assumed uncorrelated between pixels, = crfSij, 
for pixels i and j, and noise RMS equals to er^. 

Let uj denote a set of cosmological parameters. Our 
goal is to compute the full joint posterior P(w|d), which 
is given by P(u)\d) oc P(d\uj)P(uj) = C(lu)P(lu), where 
£(w) is the likelihood and P(w) a prior. For a Gaussian 
data model, the likelihood is expressed as: 



a x ,i m = J dfieVim(e) J dk<5(k) (^^j {-i) l Pi(k-e)@ ( x}(k). 

(1) 

Here, Y/ m (e) denotes the spherical harmonics, P(k) are 
the Legendre polynomials, and Q xl (k) is the Ith mo- 
ment of the transfer function of scalar modes, for either 
temperature or polarization. Further, <5(k) is a random 
variable that characterizes the initial amplitude of the 



£{u>) oc 



(7) 



where C = S + N is the total covariance matrix. 

3.1. The Gibbs sampler 

The problem of extracting the cosmological signal s 
and lo fro m the full signal by Gibbs sampling was ad- 
dressed bv l Jewell et"aTl (|2004l) . IWandelt et ail (|2004f) and 
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lEriksen et all (|2004bf ). The CMB Gibbs sampler is an 
exact Monte Carlo Markov chain (MCMC) method that 
assumes prior knowledge of the conditional distributions 
in order to gain knowledge of the full joint distribution. A 
significant fraction of the CMB data is completely domi- 
nated by galactic foreground, and about 20% of the data 
needs to be removed. This might sound trivial, but in 
reality it complicates processes as the spherical harmon- 
ics no longer are orthogonal. The Gibbs sampler solves 
this problem intrinsical ly, as the galaxy m ask becomes a 
part of the framework (|Groenebooml 2009). 

The main motivation for introducing the CMB Gibbs 
sampler is the drastically improvement in scaling. With 
conventional MCMC methods, one needs to sample the 
angular power spectrum, Cg — {ai m a^ m }, from the distri- 
bution P(Cf|d), which scales as 0(iVp ix ), where A p i X is 
the size of the covariance matrix. For a white noise case, 
the Gibbs sampler reduces this to O(N^). In other 
words, the Gibbs sampler enables effective sampling in 
the high-^ regime. 

3.2. Sampling scheme 

In order to sample from the full joint distribution 
P{Ci, w,s|d) using the Gibbs sampler, we must know 
the exact conditional distributions P(s\Ci , ui, d) and 
P(Ce,u)\s). The Gibbs sampler then proceeds by alter- 
nating sampling from each of these distributions: 

(C i ,u ) ) i+1 ^P(C £ ,^|s\d) (8) 
s i+1 ^P(s\(C i ,u;) i+1 ,d). (9) 
The first conditional distribution is expressed as: 

-is T S(w)- 1 s 

P(C t ,w\*,d)= ^r^r , (10) 

and is distributed according to an inverse Gamma func- 
tion with 2£ + 1 degrees of freedom. The remaining con- 
ditional distribution is 

P( S \C e ,u,d)o,e-^ s -^ T ^r 1 +™- 1 )(s-s) i (n) 

where s = N _1 d. In other words, P(s\Ce,u>,d) is 
a Gaussian distribution with mean s and covariance 
(S(w) _1 + N -1 ) -1 . Numerical methods for sampling 
from t hese distributions were discussed by IGroenebooml 
(|2009f ). and the details on how the polarization covari- 
ance matrix was numerically implemented can be found 
in Appendix [SJ 

4. RE-ANALYSIS OF 5- YEAR TEMPERATURE WMAP 
DATA 

lAckerman et al.l (|2007D and IHanson fe Lewis! (|2009f ) 
pointed out an error in the expression for the off-diagonal 
covariance matrix. The expression for the signal covari- 
ance matrix now includes a previously neglected 

factor of . For the ACW-covariance matrix that 

correlate scales with £ = £' ± 2, the only difference in 
contrib ution is (— i) ±2 = — 1, n egating the off-diagonal 
terms. IHanson fc Lewis! ()2009D claims that the ACW- 
signal direction in the 5-year WMAP data is located at 
the eclipt ic poles, and not at (l,b) = (11 0°, 30°), as pre- 
sented by IGroeneboom fc Eriksenl {2009) • In light of the 
new results, we perform a new full temperature analysis 



TABLE 1 

Summary of marginal posteriors from WMAP5 



Band 


t 


range 


Mask 


Amplitude g* 


Direction (£, b) 


Wl-4 
Vl-2 
Ql-2 


2 
2 
2 


- 400 

- 400 

- 300 


KQ85 
KQ85 
KQ85 


0.29 ± 0.031 
0.14 ±0.034 
-0.18 ±0.040 


(94°, 26°) ±4° 
(97°, 27°) ±9° 
(99°, 28°) ± 10° 



Note. — The values for g t indicate posterior mean and standard 
deviation. The ecliptic poles are located at ±(96°, 30°). 



of the WMAP data and investigate whether the neglected 
factor has any impact on the resulting posteriors. 

4.1. Data 

We consider the five-ye ar WMAP temperature sky 
maps (jHinshaw et al.ll2009f ). and analyze the Q-, V and 
W-bands (41, 61 and 94 GHz), where the W and V bands 
are assumed to be the cleanest WMAP bands in terms of 
residual foregrounds. We adopt the template-corrected, 
foreground reduced maps recommended by the WMAP 
team f or cosmolog ical analysis, and impose the KQ85 
masks (|Gold et alj|2009f ). which remove 18% of the sky. 
Point source cuts are imposed in both masks. 

We analyze the data frequency-by-frequency, and con- 
sider the combinations V1+V2, Q1+Q2 and Wl through 
W4. The noise RMS patterns and beam profiles are 
taken into account for each difference assembly map 
(DA) individually. The noise is assumed uncorre- 
cted. All data used in this analysis are available from 
LAMBDA 6 . 

4.2. Results 

The results from our analysis are presented in Ta- 
ble [TJ and the posteriors are shown in Figure [TJ The 
strongest detection is still present in the W-band, where 
.9* = 0.29 ± 0.031, corresponding to a 9er detec- 
tion. However, the correction term mentioned above 
cle arly has a significant effect on the signal described 
by IGroeneboom fc Erik sen (200l|). The direction and 
the significance of the detection is altered: For both the 
W-band and V-band analyses, the preferred direction is 
now located at (l,b) = ±(96°, 30°), very close to the 
north/south ecliptic poles. In addition, the significance 
of the signal in the W-band is increased from previously 
3.8c to about 9a, showing that the neglected correction- 
term has "forced" the signal away from its true direction 
- the north and south ecliptic poles. The probability 
that this direction is a pure coincidence is minimal, and 
the observed signal is therefore most likely a product of 
systematics. Another interesting fact is that the signal 
seems to be frequency dependent, with a stronger sig- 
nal in the W bands than in the V bands. Further, the 
Q-bands seems to exhibit a negative <?*, which suggests 
frequency dependence. 

5. ANALYSIS OF SYSTEMATIC EFFECTS 

Before the correction was introduced, we performed 
several tests on the independent WMAP 5-year DA 
bands showing that the direction is both existent and sta- 
ble in bands. The significance was also slightly increased, 
and we are able to cover the signal up to £ max = 700. We 

6 http://lambda.gsfc.nasa.gov/ 
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Fig. 1. — W and V-band posteriors for the temperature analysis, using t cu t B = 400 and the KQ85 mask. The north and south ecliptic 
poles are marked with a red ci rcle. Note how the posterior pe aks correspond with the ecliptic poles. The yellow circles indicate the direction 
from the previous analysis by IGroeneboom & Erikscn ( 20091) . 



now proceed by investigating various systematic effects 
as candidates for the observed signal. A visualization 
of some possible sources of systematic effects together 
with a realization of the ACW signal for comparison is 
presented in figure [H Asymmetric beams (upper right), 
noise RMS maps (bottom left) and the zodiacal light 
template (lower right). 

5.1. Impact of noise mis- estimation 

One of the possible candidates for generat in g the 
ACW-signal found by iGroeneboom fe Eriksenl (|2009f ) 
is noise mis-characterization. Previous work done by 
(|Groeneboom fe Eriksenl I2009T ) showed that correlated 
noise levels have little or no effect on the signal. How- 
ever, it might be possible that noise with incorrect RMS 
specifications could give rise to a signal similar as the 
ACW signal. We therefore perform one more analysis to 
t est noise sensitivity. 

IGroeneboom et all (|2009f ) discovered that the noise 
levels provided by the WMAP team were slightly off by 
about 0.5 — 1 %. While this error is small enough to not 
significantly affect most cosmological analyses, it is con- 
ceivable that incorrect noise levels could contribute to a 
signal similar to the ACW-model. 

We therefore simulate a VI map with 5 % incorrect VI 
noise, i.e., the noise is multiplied with 1.05 before it is 
added to the map. The analysis is done with the KQ85 
mask. The x 2 comes out about 6% above the expected 
value, recording that the incorrect noise is measured by 
the Gibbs sampler. However, the posteriors still show a 
zero-detection of the ACW-model, with an anisotropic 
amplitude of = 0.01 ± 0.05. This indicates that in- 



correct noise levels have little or no effect on the ACW- 
signal. 

5.2. Impact of asymmetric beams 

Another issue with the analysis of 
IGroeneboom fc Eriksenl (|2009f ) is whether the asymmet- 
ric beams of the WMAP detectors could h ave given rise 
to a s ignal similar to the ACW model. IWehus et al.l 
(|2009f ) established a full framework for simulating 
WMAP maps with asymmetric beams. An example of 
contribution from asymmetric beams on WMAP maps 
is presented in Figure [5J The authors also provided 
a set of 10 simulated maps with asymmetric beams. 
We now perform a Bayesian analysis on these maps, 
together with an analysis on isotropic simulated maps 
with symmetric beams for comparison. 

The test data are set up as such: We simulate isotropic 
test maps with the best-fit ACDM power spectrum, and 
convolve them with the standard symmetric V-beams. 
We then add V-band noise RMS to the maps, and analyze 
the test maps. We th en perform the sam e analysis on 
the V-band maps from IWehus et all (|2009l ). which were 
produced with asymmetric beams. Both analyses are 
done using multipoles Z max = 700 and 1™'° = 512, with 
a standard V-band setup and the KQ85 mask. 

The posteriors for the anisotropy amplitude are 
shown in Figure [3J with both having = —0.01 ± 0.05. 
It should be clear that asymmetric beams do not pro- 
duce effects in the CMB similar to the ACW-model, as 
the analysis show no trace of any signal detection. 

5.3. Zodiacal light 



■5 



ACW signal 



Asymmetric beams 




Fig. 2. — Various systematic effects compared with the ACW signal (upper left). Asymmetric beams (upper right), noise maps (bottom 
left) and the Zodiacal light template (lower right) are similar in shape to the ACW signal, and could therefore be thought to contribute to 
the ACW signal in the WMAP data. 




Symmetric beams 
Asymmetric beams 



Fig. 3. — An analysis of the same isotropic map convolved with 
symmetric (black) and asymmetric (red) beams. Note how both 
results are consistent with g t = 

In this paper, we have seen that the ACW signal in the 
WMAP data has shifted to the ecliptic poles. This in- 
dicates that the signal is most likely not of cosmological 
origin, as it is strongly aligned in the plane of the solar 
system. An interesting question is whether the ACW- 
signal is connected to the Zodiacal light. Zodiacal light is 
produced by sun-rays reflecting off dust particles sharing 
the same orbit as the earth, and share a similar overall 
structure as the ACW-signal. An illustration of a zodia- 
cal light template is presented in Figure^ together with 
the estimated ACW signal in the direction of the eclip- 
tic poles. We perform three analyses of realistic V-band 
simulations, where we co-add the Zodiacal light template 
to simulated, isotropic V-band simulated maps. In the 
first run, we add the template as it is, in the second and 
third analyses we multiply the template with a factor of 
10 and 100, respectively. In all of the analyses, the poste- 
riors resulted in zero-detections with g„ = 0.0±0.045 and 



no significant directions on the sky, with uniform distri- 
butions. We therefore conclude that the Zodiacal light 
does not have a significant contribution to the ACW sig- 
nal in the WMAP data. 



5.4. Analyzing alternative WMAP data 

iLiu fc LP (|2009al fbh have developed an alternative 
framework for building 1-year WMAP maps from raw 
data. The authors imposed stronger constraints on data 
selection, removing almost 20% of the time-ordered data. 
For instance, data for which the beam boresight distance 
from the planets are less than 7° are removed, corre- 
sponding to the the antennae main beam radius. The 
temperature map published by the WMAP team used 
a cut of only 1.5° (|Limon et aj 12001) . ILiu fc Lil (|2009aD 
also used a n extended KQ85 -mask which removes 28.3% 
of the sky. ILiu fc Lil (|2009bf ) claim that the pixels in the 
WMAP scan ring of a hot pixel are systematically cooled, 
where the strongest anti-correlations between tempera- 
tures of a hot pixel and its scan-ring appear at a sepa- 
ration angle of about 141°. Due to the anti-correlation 
of pixels and the strict data cuts , the tem perature power 
spectrum obtained bv ILiu fc Lil (|2009af ) is decreased on 
average by about 13%, causing the best-fit cosmological 
parameters to change considerably. 

In order to see whether the anti-correlated pixels in 
the WMAP stream could have contributed to the ACW 
signal in the WMAP data, we perform a full tempera- 
ture analysis on both t he alte rnative temperature map 
provided by ILiu fc Lil (|2009a| ) and the original 1-year 
WMAP temperature map. The map used in our analysis 
is the Vl-band. The RMS noise map fo r the alternative 
analysis is provided bv lLiu fcTiil (|2009a| ). while the maps 
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for the standard WMAP analysis were downloaded from 
the Lambda site. The VI beam is th e same in both ca ses, 
as is the extended KQ85 mask from lLiu fc Lil (|2009aD . If 
the ACW sig nal is detected in the WMAP data but not 
the data from iLiu fc Lil(l20f)9ah . it might be an indication 
that the WMAP team have included data that should 
have been left out, giving rise to a correlation structure 
similar to that of the ACW signal. 

Analyzing the maps up to £ max = 400, we find that 
both maps do contain a significant anisotropic signal, 
with g* ~ 0.15 ±0.10. This implies that the ACW-signal 
is most likely a more intrinsic part of the WMAP data, 
and not due to the possible anti-correlation of pixels. 

6. FORECASTS FOR PLANCK WITH POLARIZATION 

The Planck satellite will provide us with high- 
resolution CMB data of superior quality compared to 
previous CMB experiments. The Planck experiment 
also provides high-resolution polarization data, with an 
^-rangc up to 2500. As the Planck data are indepen- 
dent from WMAP data, it will be very interesting to see 
whether the ACW-signal is evident or not in the data. 
We therefore need to investigate some anisotropic prop- 
erties of typical Planck-data in order to know what to 
expect and not expect. 

In this section, we set up a high-l? temperature analysis 
with !™^ ff = 800 and a joint temperature and polariza- 
tion analysis with !™^ ff = 400. We then analyze the 
maps to obtain the posterior means and standard devia- 
tion. We continue by forecasting how the standard devi- 
ation of the anisotropic amplitu de posteriors should vary 
with multipoles £, as done by IGroeneboom fc Eriksenl 
((2009ft . 

6.1. Validation of the polarized sampler 

Before performing a full-scale analysis of simulated po- 
larized Planck data, we wish to validate our code. We 
therefore simulate a low-resolution N s i^ c = 32 map with 
E- mode data included. Assuming an anisotropic am- 
plitude of g* = 1.0, we perform both a brute-force and 
a metropolis-hastings analysis of a full-sky map with no 
beam nor noise. The resulting posteriors for the TT-case 
and the TT+TE+EE-case are shown in Figure HI It is 
worth to note that the posterior is more narrow when in- 
cluding polarization data, as there is more data available. 
A typical posterior of the estimated direction n together 
with the input TT+EE ACW-signal is seen in Figure [3 

6.2. Simulations 

We now consider a Planck simulation. We first sim- 
ulate a temperature-only ACW-anisotropic map with 
n sid0 = 1024, £ max = 2000 and 4 uto ff = 1024, with a pre- 
ferred direction pointing towards (9,<p) = (57°, 57°) and 
an anisotropy amplitude of = 0.1, us ing the best-fit 
5-yea r WMAP ACDM power spectrum (|Komatsu et al.l 
200]|. The map is convolved with a Gaussian beam 
corresponding to the 143 GHz Planck channel, and 
white, uniform noise is finally added. The beam FWHM 
for this frequency channel is 7.1', and the temper- 
ature noise RMS per A^dc = 1024 pixel is er T = 
12.2/LtK. The polarization noise RMS is op = 23.3/iK 
(|The Planck Collaboratfol3l2006fh 
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Fig. 4. — g+ posteriors for several analysis of noiseless simulated 
^max = 64 map, using both MCMC and brute-force calculations. 
Note how the polarization data narrows the distribution. 




Fig. 5. — Posterior from a simulated set with g» = 1.0. The 
original temperature ACW-signal in the input map can be seen 
in the background. Note how the estimated direction corresponds 
well with the posterior. 

6.3. Results 

We perform three analyses of the simulated Planck sky 
map. The first is an analysis on low-£ (l^^ s = 400) 
temperature data, the second high-£ ((™^ ff = 800) tem- 
perature data while the third is a low-£ analysis of 
TT+TE+EE polarization data. The results are shown in 
Table [2j where we reproduced the input parameters with 
typically g* = 0.11 ±0.025. Note how the standard devi- 
ation of the posterior is lower than for the WMAP case. 
This is to be expected, as higher multipoles £ contribute 
more to the anisotropic effect, but not significantly. This 
is due to the fact that the off-diagonal correlation terms 
in the covariance matrix have a lower values on smaller 
scales. 

We determine the standard deviation of the <?» pos- 
terior as a function of multipoles £ by simulating a un- 
convolved, noiseless isotropic map including polarization 
data using the best-fit ACDM power spectrum. We 
then analyze this map for various £, obtaining the pos- 
terior distribution for each run. The results are seen 
in Figure [6] Here, we see that cr(£*) is very close 
to a power law in £ , in good agreement with th e ar- 
guments given by iPullen fc Kamionkows"kl (|2007D and 
IGroeneboom fc Eriksenl (|2009D . The best-fit power law 
function is 

cK4igh;s*) = 0.0117 (^) 27 , (12) 
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TABLE 2 

Summary of marginal posteriors from simulated Planck data 



Simulated data 


Input amplitude 


I range 


Mask 


Estimated g* 


Low-£ TT 


0.10 


2-400 


KQ85 


0.11 ±0.025 


High-£ TT 


0.10 


2 - 800 


KQ85 


0.11 ±0.020 


Low -I TT+TE+EE 


0.10 


2-400 


KQ85 


0.10 ±0.020 



Note. — The values for indicate the posterior mean and standard deviation. 
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Fig. 6. — Estimated uncertainty in g* as a function of £ (black 
dots) and a best-fit power law function (red line) for cosmic vari- 
ance limited data. 

and this can be used to produce rough forecasts for the 
Planck experiment including polarization. For instance, 
if both temperature and E-mode polarization data are 
available up to £ = 512, then the standard deviation of 
g* is cr(512) ~ 0.001. This is generally a factor two better 
than using temperature alone. 

7. CONCLUSIONS 

We have generalized a previously developed Baycsian 
framework to allow for exact analysis of any general 
anisotropic universe models that predicts a sparse sig- 
nal harmonic space covariance matrix, including polar- 
ization data. This generalization involved incorporation 
of a sparse matrix library into the existing Gibbs sam- 
pling code called "Commander" . We implemented sup- 
port for this model in our codes, before demonstrating 
and validating the new tools with appropriate simula- 
tions including polarization data. First, we compared 
the results from the Gibbs sampler with brute-force like- 
lihood evaluations, and then verified that the input pa- 
rameters were faithfully reproduced in realistic WMAP 
simulations. 

We then considered a s pecial case of an i sotrop ic uni- 
verse models, namely the lAckerman et al.l ([20071 ) model 
which generalizes the primordial power spectrum P(k) to 



include a dependence on direction, -P(k). The equations 
we re however not complete , and t he analysis performed 
by iGroeneboom fc Eriksenl (|2009T ) has been re-done in- 
cluding the previously neglected -term. 

We then analyzed the five-year WMAP temperature 
sky maps, and presented the updated WMAP posteri- 
ors of the ACW model. The results f rom this analysis 
are in accordance with the results from lHanson fc Lewii 
(|2009( ). showing that the preferred direction is now lo- 
cated at the ecliptic poles. This suggests that the signal 
is most likely not of cosmological origin, and its origin 
must be either from within the solar system or system- 
atics. 

We have investigated four cases of systematic effects 
that share similar structures with the ACW signal. We 
have shown that neither asymmetric beams, the Zodiacal 
light, noise RMS mis-estimation nor possible pixel anti- 
correlations in the WMAP data could have given rise to 
the observed signal. 

To summarize, we have shown that there exist a strong 
anisotropic signal corresponding to the ACW signal in 
all the WMAP data that is aligned with the north and 
south ecliptic poles. The probability that the axis should 
correspond so closely to the ecliptic poles is very low, 
indicating that the signal is due to a systematic effect. 
The signal makes up more than 5% of the total power 
of the temperature fluctuations in the CMB. We have 
excluded some of the possible candidates as source of the 
ACW signal. Determining the nature of the systematic 
effect will be of vital importance, as it might affect other 
cosmological conclusions from the WMAP experiment, 
and the upcoming Planck data will clearly be invaluable 
for understanding the nature of this feature. 
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APPENDIX 
THE COVARIANCE MATRIX 

Even though we do not employ B-mode polarization data in the analysis performed in this paper, the numerical 
framework still supports B-mode polarization. In this section, we therefore describe the full TT+EE+BB covariance 
matrix including correlations. In the previous analysis, only temperature-temperature anisotropic correlations were 
considered. We now extend the framework to include polarization, such that the Fourier coefficients become: 

{ TT EE BB \ { a \ 

aim = [a lm , a tm , a lm ). (Al) 



s 



The TT+EE covariance matrix for lmax-8 





Fig. 7. — The ACW TT-EE covariance matrix (left) in the representation of Equation IA2I Diagonalizing this matrix turns out to be 
a nearly impo ssible task, forcing us to use another representation. The ACW TT-EE covariance matrix (right) in the representation of 
Equation IA3I Diagonalizing this matrix is similar to diagonalizing the TT-only ACW covariance matrix, and is more efficient. 



The covariance matrix C, 



tm.t' 'm' 



can be expressed as: 



/TT TE TB\ 
Clm,£'m> = TE EE EB 
\TB EB BB) 



(A2) 



The existing framework for sampling anisotropic universe models in FORTRAN was then altered to allow for po- 
larization data, and whether polarization is used is flagged through a parameter file. The off-diagonal TT+EE+BB 
anisotropic covariance matrix is presented in figure [7] Note that the BB component is zero in this plot. However, this 
straight-forward representation of the full covariance matrix is too naive: performing a Cholcsky-factorization (diag- 
onalizing) of this matrix for high £s arc nearly impossible. Diagonalizing a matrix is more efficient when off-diagonal 
elements are close to the diagonal. However, the (TT,EE,BB) representation of the matrix in figure [7] gives rise to 
elements spread around the full matrix. Typically, Cholesky-factorization for such a TT-EE-BB representation breaks 
down for Z max = 64 due to the dense structure of the upper-triangular decomposed L-matrix. 

To overcome this problem, we ope rate with a different representation of the (TT,EE,BB)-matrix. Instead of building 
the matrix as presented in equation IA2| we choose a different way of expressing the matrix: 



(A3) 



with corresponding ae m s 





1 TTqq TEqq 


TB o 


TT 10 .. 


EE n0 \ 




TEoo EEoo 


EB o 


TE W . . 


EB n0 




TBqq EB o 


BB aa 


TB W .. 


BB n0 




TT m TEqi 


TBoi 


TT U .. 


EE n i 




\TB 0n EB 0n 


BB an 


TBi n . . 


BB nn J 




( a 00i a 00i a CK)i 


T 






a-im — 


a 01 , . . . 





(A4) 



As the EE and BB correlations share the same structure as the stand-alone TT, the complete covariance matrix will 
in this representation resemble the original three-banded covariance matrix. The elements are now much closer to the 
diagonal, solving the problem of inefficient diagonalizing. The matrix representation is depicted in figure Note that 
this representation is only used when multiplying the matrices with vec tors and performing Cholesky decompositions. 
Withing the rest of the framework, the a^ m s are treated as in equation lAll 
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